Collective motion due to individual escape and pursuit response 
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Recent studies suggest that non-cooperative behavior such as cannibalism may also be a driving 
mechanism of collective motion. Motivated by these novel results we introduce a simple model of 
Brownian particles interacting by pursuit and escape interactions. We show the onset of collective 
motion due to escape and pursuit response of individuals and demonstrate how experimentally 
accessible macroscopic observables depend strongly on the ratio of the escape and pursuit strength. 
We analyze the different impact of the escape and pursuit response on the motion statistics and 
determine the scaling of the migration speed with model parameters. 

PACS numbers: Valid PACS appear here 



The emergence of collective motion of living organisms, 
such as exhibited by flocks of birds, bacterial colonics 
or insect swarms is an ubiquitous and fascinating self- 
organization phenomenon in nature, which still bears 
many open questions. A common explanation for the 
emergence of collective motion in a wide range of ani- 
mals is that it serves as a protection mechanisms against 
predators. Recent experimental results suggest a novel 
mechanism driven by cannibalism which may, surpris- 
ingly, facilitate collective motion in mass migrating in- 
sects [l|, [IJ. Insect swarms can extend over many kilo- 
meters and can have a devastating impact. Locusts, for 
example, can invade up to one fifth of the Earth's land 
surface and are estimated to affect the livelihood of one 
in ten people in the planet Q . 

The phenomenon of swarming in general has attracted 
scientists from a wide range of disciplines with different 
scopes and perspectives [|| . In recent years it also be be- 
came the focus of an increasing number of publications in 
the field of statistical physics, non- linear dynamics and 
pattern formation. These contributions enhanced signifi- 
cantly our understanding of collective motion in systems 
of self-propelled particles (SPP) by discovering univer- 
sal scaling laws and phase-transition like behavior, and 
offered new stimuli to the theory of non-equilibrium sys- 
tems 0, i, 0, 1, i, M, El HUSH- Recent examples 
of the ongoing research are a special journal issue dedi- 
cated exclusively to active motion and swarming (l5j or 
the work by Grossman et. al. [l6| . where the authors 
demonstrate the emergence of collective motion of SPP 
interacting via inelastic collisions. 

In this letter we investigate a simple but generic model 
of individuals with escape and pursuit behavior which 
may be associated with cannibalism. We show the onset 
of collective motion in the absence of explicit velocity- 
alignment interaction, which to our knowledge is an inte- 
gral ingredient of most related models. Our work demon- 
strates a strong dependence of macroscopic dynamics on 



the relative strength of individual escape and pursuit re- 
sponse. Directed translational motion in our model is 
a strictly collective (but not cooperative) behavior and 
may be therefore termed group-propulsion. Our model 
offers a novel perspective on possible mechanisms of on- 
set and persistence of collective motion and the resulting 
migration patterns in nature and represents an interest- 
ing example of pattern formation and phase-transitions 
in non-equilibrium systems. 

We model an individual organism as an active Brow- 
nian particle in two dimensions (d = 2) with an internal 
energy depot (see [13] for details). This additional degree 
of freedom describes the energy budget of individuals de- 
termined by their uptake of nutrients, internal dissipation 
to maintain body processes and conversion of energy into 
energy of motion. It allows individuals in our model to 
increase their speed in reaction to external stimuli by 
conversion of internal energy into energy of motion. For 
simplicity we assume throughout this work that at all 
times there is a surplus of internal energy which allows 
us to neglect the explicit treatment of the energy balance 
and focus on the spatial dynamics only. 

Each individual (particle) obeys the following 
Langevin dynamics: 



Ff + s/lD^ii (1) 
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The first term on the left hand side of the velocity equa- 
tion (fT]) is a friction term with friction coefficient 7 and an 
arbitrary power dependence on velocity represented by 
a = 1, 2, 3, The response of individual i to other in- 
dividuals is described by an effective social force F| . The 
last term is a non-correlated Gaussian random force with 
intensity D v . A solitary individual (Ff = 0) explores its 
environment by a continuous random walk, where the 
individual velocity statistics are determined by 7, a and 
D v . The parameters are given in arbitrary time and space 
units T and X: [7] = X x ~ a T a - 2 [D v ] = X 2 T~ 3 . 

The finite-size of individuals is taken into account by 
fully elastic hardcore collisions with a particle radius Rhc 
(for details see f\l§). 

Motivated by experimental observations [H, we intro- 
duce the following response mechanisms: If approached 
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from behind by another individual j the focal individual 
i increases its velocity away from it in order to prevent 
being attacked from behind. We refer to this behavior 
as escape (e). If the focal individual "sees" another indi- 
vidual up-front moving away, it increases its velocity in 
the direction of the escaping individual. We refer to this 
behavior as pursuit (p). No response in all other cases. 
The response of an individual is determined the follow- 
ing decision algorithm: A) Is there another individual 
within my sensory range l s > Rhc] B) If yes, is it in front 
or behind me, and C) does it come closer or does it move 
away. 

Based on the above considerations we write Ff as a 
sum of an effective escape and an effective pursuit force: 
Ff = ff + ff with 

{ >= W e E Av i<*& - »*)0(-Vif jiM-vj,?*) (2a) 
e j 

ff = |FE Av i«*& ~ rumviijiMvjitji) (2b) 
p j 

where Xp,e > are the corresponding interaction 
strengths, Avji = (v^r.,-,-)?^ is the relative velocity of 
particle j with respect to particle i, with Vjt = Vj — Vj and 
*ji = ( r i ~ r i)/\ r j ~ r i\- The Heavisidc functions 6 reflect 
the conditions for the escape and the pursuit response. 
Both forces are normalized by the respective number of 
individuals which the i-th individual responses to: N e 
is the number of individuals which fulfill the escape re- 
sponse conditions and N p the corresponding number for 
the pursuit response. 

The symmetry of the introduced interaction is broken 
in several ways: the interaction acts only on one of the 
interacting particles (action^ reaction); the interactions 
are direction selective - the particles distinguish between 
their front (vj • Tji > 0) and their back (vj • r^, < 0) and 
between approach (Vjj-r^ < 0) and escape (v 3 -j- r ji > 0); 
the strength of interaction to the front and back may be 
different (xe 7*= Xp)- The most important property of the 
interactions is their anti-dissipative nature with respect 
to kinetic energy. Note that Ff leads only to accelera- 
tion of individuals and is analogous to the autocatalytic 
machanism proposed in Bazazi et al. Q ■ 

Throughout this letter we will discuss our numerical re- 
sults in terms of the rescaled density p s = Nl^/L 2 , where 
N is the total particle number, l s the interaction range 
and L the size of the simulation domain. All simulation 
results were obtained with periodic boundary condition. 
We will restrict here to the case of moderate noise in- 
tensity D v < 1 and focus on the system behavior with 
changing density p s , which may easily be controlled in ex- 
periments. Although the system shows interesting global 
behavior with increasing noise D v , related to phase tran- 
sitions known from systems of SPP [!, H, [H| , the corre- 
sponding analysis is beyond the scope of this letter and 
will be discussed in a forthcoming publication. 

Numerical simulation reveal that irrespective of the de- 
tailed model parameters, the pursuit and escape inter- 
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FIG. 1: Typical spatial configurations and particle velocities 
(small arrows) for pursuit only (p), pursuit+escape (p+e) and 
escape only (e) cases at different particle densities p s = 0.30, 
1.25, 2.25. The direction of the large arrows indicates the 
mean migration direction and their length the migration 
speed U . For the escape only case at low densities the mean 
migration speed vanishes U ~ 0. 



actions lead to global collective motion at high particle 
densities p s and moderate noise intensities D v (Fig. [1]). 
At low p s however we observe a very different behavior 
in dependence on the microscopic details of the model, 
where the velocity statistics and spatial migration pat- 
terns depend strongly on the relative strength of the es- 
cape and pursuit interaction \y an d Xe- For \ e > 
and Xp ~~ * with increasing p s a transition between a 
disordered state, with vanishing mean migration speed 
(U) = I E, V i\/N = and an ordered state with (U) > 
takes place. This resembles similar transitions reported 
for SPP with velocity alignment (Fig. UJ^ ||. With in- 
creasing Xp the transition shifts to lower p s until it van- 
ishes. For Xp > and x e — > there is no dependence of 
(U) on p s . 

In order to understand the dynamics we investigate 
the influence of escape and pursuit interactions indepen- 
dently, by analyzing the extreme cases: Xp — 0> Xe > 
(only escape) and Xe = 0, Xp > (only pursuit). 

In the escape only case the particles try to keep their 
distance with respect to individuals approaching from be- 
hind. To the front only interactions via the hardcore col- 
lisions take place. At low p s after an escape response the 
probability of interaction within the characteristic time of 
velocity relaxation vanishes and the particles are able to 
reorient themselves (disordered state). As p s increases 
the frequency of escape interactions increases and the 
particles are able to correlate their velocities on several 
interaction length scales but subensembles may move in 
different directions. We observe a transition to the or- 
dered state via an active fluid like state (p s « 1.25; Fig. 
[2^). In the ordered state all particles are able to correlate 



3 




00" x*x' 
X x 


p s =0.24 


(c) go 

Pa 
/V>Dqd 
/ « D 8° 

/ Br 


p s = 1.68 



FIG. 2: (a) Mean global velocity (U) for escape-only (O) 
Xe = 10, Xp = 0, pursuit-only (x) \e = 0, \v = 10 and 
symmetric escape+pursuit (□) Xe = Xp = 10 over particle 
density p s obtained from numerical simulations with periodic 
boundary conditions (7 = 1, D v = 0.05, a = 3, Rhc = 1, 
l 3 =4). Only pure translational solutions were considered 
and the errorbars represent one std. deviation; Particle speed 
distribution P(v) for the different interaction types in compar- 
ison with the analytical solution for non-interacting Brownian 
particles (solid line) at p s = 0.24 (b) and p s — 1.68 (c). 
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FIG. 3: (a) Snapshot of a single large vortex formed from ran- 
dom initial conditions for the pursuit only case. The arrow 
indicates the rotation direction; Comparison of numerically 
obtained average pair velocities (it) for a = 1 (O) and a = 3 
(□) with the result of Eq. U (solid lines): (b) (it) over fric- 
tion coefficient 7, (c) noise intensity D v , (d) and interaction 
strength x- ( e ) { u ) vs - friction function exponent a. Here 
we distinguish two cases A > 1 (O) and A < 1 (□), with 
A = y/D vX /n/j. 



their direction of motion. At all p s we obtain spatially ho- 
mogeneous distribution of particles. The transition-like 
behavior is also reflected in the particle speed distribu- 
tion P(v). At low densities P(v) corresponds to the an- 
alytical result obtained for non-interacting particles (Fig 
Wp), whereas at high p s the maximum of the distribution 
shifts to higher speeds indicating a transition from pure 
random walk to directed translational motion (Fig[2j:). 

In the case of pursuit-only interaction the dynam- 
ics change dramatically. Already at low p s we observe 
a highly inhomogeneous state: initiated by formation 



of small compact particle clusters performing coherent 
translational motion. As there is no escape interaction 
the density of the clusters is only limited by the hard- 
core radius. At moderate noise intensities the clusters 
are highly stable and a process of cluster fusion can be 
observed where larger clusters absorb smaller clusters 
and solitary particles. The dominant stationary con- 
figuration with periodic boundary condition, and mod- 
erate noise, is a single large cluster performing trans- 
lational motion (Fig. [T|). The migration speed (U) in 
Fig. [5^ is given by the mean speed of a single cluster 
( u ) = I Sieciustor v * I /^cluster, which for large clusters be- 
comes independent of the cluster size and therefore inde- 
pendent of p s . The same holds for P(v) as shown in Fig. 
[5j3,c. An intriguing feature of pursuit-only is the possi- 
bility of the formation of large scale vortices out of ran- 
dom initial conditions due to collisions of clusters moving 
in opposite directions. After a nucleation a vortex may 
grow by absorbing smaller clusters leading to a single ro- 
tating structure (Fig. (3£i) with life times exceeding 10 3 
time units. Preliminary results on vortex-stability indi- 
cate a monotonous increase of stability (i.e. life time) 
with size (not shown). The emergence of vortices in our 
model is in particur remarkable because so far they have 
only been reported for systems of SPP with confinement, 
or attracting potential, respectively EH, H(| • Here the 
pursuit interaction acts in a sense as both: a propulsion 
mechanism and an asymmetric attraction. 

The analysis of the dynamics shows that both interac- 
tions — escape and pursuit — lead to collective motion 
but have an opposite impact on the density distribution. 
Whereas the escape interaction leads in general to a ho- 
mogenization of density within the system, the pursuit 
interaction facilitates the formation of density inhomo- 
geneities (clusters). This leads us to the insight that the 
actual escape+pursuit dynamics where XpiXe > is a 
competition of the two opposite effects with respect on 
the impact on the particle density. The stability of mov- 
ing clusters in this simple model is determined by the 
relative ratio of the interaction strengths. In general for 
the escape+pursuit case at low p s we observe fast for- 
mation of actively moving particle clusters with complex 
behavior: fusion and break up of clusters due to clus- 
ter collisions as well as spontaneous break up of clusters 
due to fluctuations. The weak dependence of the par- 
ticle speed distribution P(v) on p s combined with the 
clear deviation from the non-interacting case at low p s 
(Fig. [2t>, c ) shows that the increase of (U) with p s for 
escape+pursuit originates in an alignment of individual 
cluster velocities. 

In order to determine the scaling of (U) with model pa- 
rameters, we consider the smallest cluster which shows 
directed translational motion: a particle pair (1,2). 
We assume particle 1 being in the front of particle 2, 
Xe = Xp = X an d |ria I < h a t all times. Through 
a transformation of Eq. [T] into polar coordinates with 
Vj = (vi cos fi, Vi shxipi), where (fi is defined as the 
angle between Vj and fi2, it can be shown that for 
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— 7r/2 < ipi < ir/2 (i = 1, 2) the escape and pursuit inter- 
action lead to an increase of either V\ or w 2 in order to 
harmonize the speed of the slower particle with the faster 
one. The acceleration is counterbalanced by the frictional 
force and results in a non- vanishing translational velocity 
of the particle pair. In addition the interaction stabilizes 
the translational motion along fi2, i.e. (ipi) — > 0. After 
the system relaxes to a stationary state (f 12 varies slowly 
in time) we end up with effectively one-dimensional trans- 
lational motion of the particle pair. 

The evolution of the mean speed of a particle pair in 
this one-dimensional situation uid = (i>i + t>2)/2 with 
vi ks vi can be approximated as 

^^-^ d +\x{\Svi\) ld . (3) 

The second term on the right hand side of Eq. [3] ac- 
counts for the acceleration of the particle pair due to the 
escape+pursuit interaction with Svi = u — «j. The factor 
1/2 accounts for the fact that at a given time only one 
of the particles accelerates. 

The deviations of individual particle speed from the 
mean speed result from the action of the non-correlated 
stochastic forces. We approximate the expectation value 
of the speed deviations (\8v\)id, by considering the speed 
deviations as discrete increments taken from a Gaussian 
distribution with zero mean and variance af d = 2D v t 
(Wiener process). Replacing r by the relaxation time of 
the interaction x _1 yields: = 2yjD v /-KX- The 

stationary velocity of a particle pair can be calculated 
from © to: 

This result is in excellent agreement with numerical sim- 



ulations of individual particle pairs (Fig. [3]) for wide pa- 
rameter ranges. The scaling in Eq. 2] agrees with the 
measurements of the average speed of large clusters at 
moderate D v and with results of a mean field approxi- 
mation which will be discussed elsewhere. 

In summary, we have presented an individual based 
model for the kinematic description of large groups of in- 
dividuals, where each individual responses to others in its 
local neighborhood by escape and pursuit behavior. The 
response itself is described by an effective social force 
acting on each individual and is motivated by recent ex- 
perimental results on mass migrating insects. 

We have shown the onset of collective motion due to 
the escape and pursuit interaction. The analysis of the 
model dynamics shows that the macroscopic behavior, 
which can be observed in experiments, such as migration 
speed vs. density or the spatial migration patterns de- 
pend strongly on the relative strength of the escape and 
the pursuit behavior of individuals. Furthermore we were 
able to obtain the right scaling of the migration speed 
with model parameters which is confirmed by numerical 
simulations. 

On the one hand, recent experiments on marching in- 
sects suggest that escape dominates their marching be- 
havior [2j . In this case our model predicts a phase tran- 
sition like behavior of the mean migration speed (U) in 
dependence on the density p, which is supported by previ- 
ous results [2l|. On the other hand, the coherent moving 
clusters and vortex structures observed in our model for 
the pursuit only case resemble observations of fish schools 
[22LI23T] and suggest the relevance of our model to a wide 
range of swarming phenomena in nature. 
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